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Synopsis 


«  •  ** 

The  partial  penetration  of  a  pumping  v/ell  is  important,  not  only  because 
it  affects  the  efficiency  of  the  well  buV.also  because  it  influences  the 
distribution  of _  head  in  its  vicinity,  modifying  the  drawdowns  observed  in 
nearby  observation' wells .  In  running  pumping  bests  to  determine  trans- 
missibility,  the  observed  drawdowns  must  be  adjusted  for  this  factor*  This 
article  outlines  the  method  of  making ‘the  necessary  adjustments.  Graphs  and 


a  table  are  given  for  partial  penetration  in  homogeneous  artesian  beds, 
supplementing  the  cursory  treatment  given  this  problem’  in  *'ate r^Supply  Paper,  ..t*  *-■ 
887. 
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The  partial  penetration  of  a  pumping  veil  is  a  matter  of:igr eat  practical 
interest.'  This  is  true  first  because  the  productivity  of  a  well  depends 
largely  on  the  degree  to  which'  it  penetrates  a  water-producing  bed.  .Second, 
t  the  distribution  of  head  in  the  -  vicinity  of  a  partially  penetrating  well  -is 
modified  according  to  the  degree  of  penetration,  and  therefore  water  levels 
'  observed  in  nearby  observation  Wells  during  pumping  tests  should  be  adjusted 
for  that  effect  before  an  attempt  is  made  to  determine  the  t ransmissibility 
or  .permeability.  \  ,  /  : 


The * problem  of  a  partially  penetrating  well  was  analyzed  in  great  detail 
by  Mu  skat  in  1932,  (See  Physics, '  vol.  '2,  pp-.  329-364;  also  "The  flow  of 
homogeneous  fluids  through  porous  media," ; McGraw-Hill/  1937,  ppv  263-286.) 
he  gave  .an  .equation,  for.  the  discharge  of  a  partially  penetrating  well  in 
terms,  of  sand-thickness,  permeability,  fractional  penetration,  and  total 
potential-drop .  f  r'om  the  assumed  effective  external'  radius  of  the  system  to 
the  known  internal  radius.  (See  equation  6,  p.  274-' of  Muskat’s  book.)  In 
1933  Kozeny  (tlasserkraft' und  Wasserwirtschaft-,'  vol*:  28,-  p.  101)  published  a 
.somewhat . simpler  empirical  expression  Summarizing  Muskat’s  analysis  to  a 
sufficient  approximation  formsny  purposes.  In  our  notation  Kozeny’s 
empirical  formula  is  ; 


1/2 


Q  =  2rrTswO<  Cl  +  7(rw/2Grm)  cos(lf  0(/2)J  /loge(re/rw) 


(D 


or 


1/2 


Q.  -  Q  9<  [x  +'7(r /  cos(rro:/2)  =  Q  c*/c 


(2) 


Here.Q  is 'the  discharge  of  a’V/ell  of  radius  r  penetrating  a  fraction 
of  cj  bed'  of  thickness  m  and  tranmissibility  f;  ^  is  the  discharge  of  a 
well  of  the  -same  'effective  radius ‘penetrating  the  bed  •  completely.  ’The  draw¬ 
down  -iri  the  'well  (i .e at  n.  distance  r  .  )  •  is  s  ,  The  drsmvdov.'n  at  the 
vi  rtu'al -external  radius  of  the  system,  re,  is  Yero.  The  quantity  in  "bracket® 
.^/abbreviated  as  (l/c).  ■  •  » 


and 


,  The  apparent  transmissibility,  referred  to  the  total  thloknees  of  themed 

determined  from  observations /of-  the  drawdown  m  a.  pumping  v  ,  S 

v  *  v. 

(3) 


T* 


Similarly  the  true  transmis sibility  is 
■T  »  Q  log  (r  /r  )  /21Y s 

,  t/t«  =■  q0/q  3  c/oc 


w 


(4) 

(s) 


The'  corresponding  ratio  of  t roe  to  apparent  permeability,  the  latter 
being  referred^however,  to  ' anequivalent  bed  of  thickness  «»  instead  of  to 
the  total  thickness,  is 


C<Q  /Q  =  C 


(6) 


This  "correction  factor"  is  given  by  Wenzel  (Water-Supply Paper  887,  p. 
,on  i  •  i?4')  It  should  be  emphasized  that  this  correction  or 

109,  equation  12  K  ,  ,  permeability  determined  from  .the  drawdovm 

adjustment  is  only  for  values  of  equation  6  (or  Wenzel's  equation 

observed  inside  a  pump_ing  weU.  Moreover,  equation  ^  K(vh„aiiv 

124)  does' not  lead  to  a  soiutTon  of  the  problem,  as  is  unKnown. 

'  .  ,  ,  • j  _ t.T,-,  o  Va+ *  c  pmiption  5.4(6),  o] 


124^  does  not  lead  to  a  solution  oi  xne  prourwu,  t7  -  fr^t^ined 

i  may  be  determined  from  Muskat’s  equation  5.4(6),  or  it  may  be 

.  .  r\  _ _ 


from  equation  2  above. 


r'  7‘h \rn “ 

Muskat,  p.  2bo.;  it  snouiu  ut=  F  stratified  and  therefore 

anisotropic  as  regards  permeability- 


•'  The  critical  distance  referred  to  is  to  be  multiplied  V  the  ,f  # 

root  of  the  ratio  of  horizontal  to  verxical  ri°ontal  in  the 

aed  is  on  the  average  sixteen  times  as  wfond  a  distance 

vertical  direction,  then  pure  y  r  accordingly  care  should  be  taken 

equal  to  about  8  times  the  mining  the  trans- 

to  adjust  drawdowns  observed  within^  that  (  •  .  _  - 

missibility  in  the  case  cited. 


In  order  to  make  the  adjustment  for  P®r^ial  ^et[^°fore going  relations 

so  know  the  distribution  of  head  throi'g  ou  ^  the’head-dif  fere  nee  between 
ire  inappropriate  because  they  relate  n.re  y  to^he ad.  ^  &  p&rtially 

the  inflow  and  outflow  surfaces.)  e  *  equation  5.3(9)  of  Muskat’s 

penetrating  veil:  is  given  byequaion  .  Jpot  COJnpUtations  for  small 

book.  The  first  mentioned  equation  con  -  ti  5.3(9)  contains  two 

values  of  rj  the  :'se,ond  for  large  jaiues  of  r.  Ration 
terms,  one  for  the  logarithmic  distribution  of  head  for P^ly 
the  other  for  the  difference  between  the  attention 

logarithmic  dist riDution.  tor  convenienc  ‘  .  d  along  the  bottom  of 

to  the  distribution  of  head  along  the  top  of  ^  f d  of  ®bs6rYation  well, 

the  bed,  such 'as  given  by  drawdown  readings  in  one  series 
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or  piezometers  just  reaching  to  the  top  of  the  bed  and  another  reaching 
to  the  bottom  of  the  bed  and  open  only  at  that  level.  The  divergence  of  the 
head  from  a  pure  logarithmic  distribution  at  a  distance  r  from  the  well  is 
given  by  ,  .  ....  m 

As/Cq/SiTT)  =  (2/rfoO^-^  ^(^fl)nKo(n-{t r/m)sin(nfi'c*)^]  /n  (v) 

the  plus  sign  being  taken  for  the  drawdown  distribution  along  the  top  of  the 
bed,  the  minu?  sign  for  the  bottom  of  the  bed. 

'  ’  *■  •  d  -  -  *  ;*•*  ••  *!  v  i 

In  equation  7  the  K  stand’s  for  the  modified  Bessel  function  of  the 
second  kind  of  zero  order.  The  other  symbols  have' the  same  meanings  as. were 

assigned ,  previously.  '  /  ....  .  ..  .  v.  ‘  " 

•  ..»>.-**  •  « 

Figures  1  and  2  were  computed  from  equation  7.  Figure  I,,  gives  the 
drawdown  adjustment  for  the  top  of  the  bed,  that  is,  the  difference  between  the 
drawdown  at  the  top  of  the  bed,  at  a  given  distance  for  different  values,  of 
fractional  penetration,-  and  the.  drawdown  at  the  same  point  with  complete, 
penetration  of.  a  well- discharging  at  the  same  rate.  Similarly,  .  figure  2  gives 
the  draw-down  adjustment"  for  the  bottom  of  the  bed,  or  the  difference:  between 
the  drawdown '.at  the  bottom  of  the  bed,  at  a  given  distance  for  different 
values  of  fractional  penetration,  and  the  drawdown  at  the  same,  .point  with 
complete  penetration  of  the  pumping  well,  the  discharge  again' remaining  fixed. 

:-aC  J.  4  '  '  ‘  . 

It. is  assumed  that  the  well  penetrates  the  bed  partially  from  the  top. 

In  that  case,  because  of  the  convergence  of  the  stream-lines  upon  th©:  screen, 
the  drawdown  will  obviously  b6  greatest  along  ' the.  top  of  the  beci  and  least 
along  the  bottom.  If  the’ well  is  screened  only  at  the  bottom_of  /the  bed  the 
problem  is  similar  though  inverted,  assuming  of  .course’ ‘that  no  part  of  the 
bed  is  unwatered.  If  the  screen  is  in  some  intermediate  position  the  curves 
in  figures  1  and  2  do  not  apply.  Nor  do  they/ apply'/  'strictly  speaking,  to 
unconfined  beds  7 •  ;  j  - . 

An  example  may  serve  to  illustrate  the  application  of'  the  drawdown 
adjustments.  Assume  that*  a  'well  with  an  effective-  r&dius  of  1.0  foot  and 
with  a  screen/' set  in  the  'tOp  40  feet  of  an  artesian  bed  100  feet  thick  has 

■.  t  -  .  j  ,  *  •  • 

been  pumping  at  a/ -constant  rate  of  840,000  g.p.d.  'long  enough  to  establish 
steady  radial  flow;  beyond ‘the  most  distant  of  three  observation  wells.  The 
observation  we  11s  are  respectively  50,  100,  and  150  feet  from  the  pumping 
Y(§11.  _  They  just -.tap  the  top  of  the  bed,  as  shown  in  figure  4.  The  drawdowns 
observed  in  those  wells  are  5.02  feet,  3’.18  feet,  and  2.30  feet  respectively,  . 
as  indicated  by  the  open  circles  plotted  in  figures  3  and  4.  The  drawdown 
in  the  pumping  veil  is  18.4  feet.  The  screen-loss  is  negligible. 


■Assume  that  the  bed  is  homogeneous  and  isotropic  as  regards  permeability. 
The  problem  is  to  determine  the  transmissibility  from  the  observed  drawdowns. 
The  procedure  is  as  follows: 


First,  plot  the  drawdowns  against  the  logarithm  of  the  distance,  as  on 
figure  3.  Draw  a  straight  line  (line  l)  as  closely  as  possible  through  the 
points  plotted.  Make  a  preliminary  determination  of  transmissibility,  uping., 
the  relation 


T  =  2.30  Q/2-rrAs. 


.v.’.v  ’ 


(8) 


«-  V  3 
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••.where  'As  is  the  drawdown  difference  over  one  log  cycle,  v  •  • 

•  ■  ■  •  •  ;  'f  ■  •  • 

Substituting,  one  gets  -  f 

T  =  2.30  x  840,000  gpd/2 TT  (9 . 1-r*  3.2)  ft  =  52,000  gpd/ft, 
and  ;  :  :  ‘  ! "  C)  '*  <,  ? 

Q/2frT  -  As/2.30  =  5i 9  ft/2*30'*  2.56  ft. 

r  -  \  .  •  -  :  ;  ,’;/J  . 

Assuming,  as  a  first  approximation >  t.hat  the  above  value^.of  T  ife  correct,' 

make  the  following  adjustments,  remembering  that  o(  =  0.40. 

•  . ,  ....  ...  n:; 

;  Tie  1 1  , 

Distance,  r,  ft 
Tt  r/m 

Observed  a2.  .  ft 

^  O  Tl  y  | 

■•Adjustment,  from  figure  1 

^stop’  ^  v  : 

Ad ju sled  d rawd own ,  ft 

•  -.  Plot  the  new  values  of  adjusted  drawdown  (x*  s  on  figure  3).  Draw  a 

.straight  line  through  the  points  (line  II ).  Make  a  second- trial  deter- 
’-mination  of  T-:  -  .•••• 


.  1 

•  2  Cm 

.3  • 

50 

.  ,J  •  ■  100 

•..-150  ■ 

1.57 

"3.14 

4.71 

5,02 

.3.18 .-. : 

2.30 

0.31 

0.0,45 

0.008 

0.80 

0.12  •;  , 

0.02  ■ 

4.22 

3.06  ; 

-  2,28.  , 

r-and 


T  =  2.30  x  840,000  gpd/2rf(7.1  -  3.l)  ft  =  77,000  gpd/ft. 


V  Q/2ifT"=  As/2.30  =  4.0  ft/2.30  «  1.74  ft.  ,  w 

*  :v  Readjust  the  observed  drawdowns,  proceeding  as^before:  rr 


ryjell  -  ■  -  ■ 

Observed’’ S£0p}  ft 
Adjustment,  from  figure  1 
/^stop,  ft 

: Adjusted  draV/down,  ft 


1  .  -2  ;3 

5.02  .  3 A 18  >  -.2.30 

0.31  Q.045  .  0.008. 

0.54  0.08  0.01 

4.48  ..  3.10  2.29 

■•u/'"  a;. 


1  H  •  •  Again  plot  the  new  values  of  a  d justed,  drawdown  (+1  s  on  figure  3).  Draw 
a  straight  line,  through  the  points,  (line  III).  Make  a .  third. .trial  deter- 
'  mindtibn:  ’ 


and 


•  T;*.2i30  x  840,000  gpd/2tt(7.7  -  3.l)ft  =  67,000 -gpd/ft ,/ 

Q/2TTT  =  As/2.30.  =  4.6  ft/2.30  »  2.00  ft.  •  ...” 


Carrying  the  procedure  through  a  fourth  time  one  gets  line  IV,  from 
which  T  =  69,000  gpd/ft  and  Q^iTT  =  1.94  ft.  This  is  close  enough  to  the 
probable  true  value.  Further  computations  are  not  called  for. 

As  a  further  check,  however,  one  may  determine  the  drawdown  that  would 
occur  in  a  completely  penetrating  well  from  Kozeny’s  formula  (equation  l). 
Solving  that  equation  for  .  '  . 

(Q/2 ITT )loge (re/rw)  and  substituting  in  it  the  proper  values,;  ^one  should  get 
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•  •  y  i".'  \  .  * 


i  • : 


drawdown  that  is checked ;  "by r  the  intercept  of  line  IV  on  the  line  r  =  r  rr 
1  ft  in  figure’  3;;  "’The  hypothetical  drawdown  comes  out' to  be  (swo('/c).' 
'this  case,  -■  — . 

^7 (cos  0.628>/8.94]  »  1.634 


In 


1  + 


(l/c).=.  I"i  7  Lcos(O.En)  j  / V  80 

and  therefore.  . .  . 

(swo^/c)  =  18.4  ft  x  0.40  x  1.634  -  12.0  ft,  which. is  verified  in 
figure  3.  ,;..  V' '  ; 

The  lower  dashed  line  in' figure  3  gives  the  probable  distribution  of 
drawdovm  ?.long  the  top  of  the  bed.  The  drawdowns  that  would  have  been 
observed  in  wells  1,  2,  and  3  had  they  reached  to  the  bottom  of  the  bed  and 
had  they  .been  open  only  at  that  level  may  be  determined  by  making  a  further 
adjustment,  using  figure  2; 

Hell 

Adjusted  -drawdown,  ft 

(from  line  IV,  figure  3) 

Adjustment,  from  figure  2 
AsbQt,  ft 


’bot  * 


ft 


1 

2 

3, 

4.42 

3.09 

2.28 

-0.28 

-0.044 

-0,008 

-0.54 

-0.09 

-0.01 

3.88 

3.00 

2,27 

:  'For  wells  1  and  2  these  points  are  plotted  as  closed  circles  on 

figure  3.  The  distinction  could  not  be  made  in  plotting  the  data  for  well 
3  and  therefore  that  closed  circle  is  omitted. 

The  upper  dashed  line  in  figure  3"  gives  the  probable  distribution  of 
drawdovm  along  the  bottom  of  the.  bed.  ,  The  drawdovm  at  the  point  of 
stagnation  underlying  the  well  may  be  determined  quite  readily  to  aid  in 

’drawing  the  drawdown-distribution  curve.  • 

•  -  S  ,  --•'i 

■4  -  /  .  r 

’  /■  .  p 

Table  1  gives  the  difference  between  the  drawdown  at  the  point  of 
stagnation  and  the  drawdovm  at  a  distance  of  4m  in  units  of  (q/2TtT)  for 
different  values  o.f  O  under  two  alternative  assumptions  regarding  conditions 
at  the  well-face.  The  data  in  the  second  column  were  computed  by  the 
writer  from  Muskat?  s  equation  5.3(8)(,  which  is  based  upon  the  assumption 
that  the  flux  thro  ugh  the  face  of  the  well,  that  is,  the  screen-velocity, 
uniform.  The  data  in  the  third  column  were  made  available  to  the  writer 
through  the  courtesy  of  Dr,  Muskat  and  the  Gulf  Research  and  Development 
Cqmpany.-  They  were  obtained  by  adjusting,  more  or  less  by  cut'  and  try, 
the  flux-density  distribution  so  as  to  give  essentially  uniform  potential 
distribution  over. the  well-face. 
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Table  1. 


\ 


Fractional 
penetration 
•  Q4 

0.05 
.  .10 
.20 
.25 
.30 
.40 
.50 
.60 
.  vo 
.7  5 
.80 

.90  . 

l.oo  *':\m 

a)  by  interpolation 

^  .  i 


Drawdown  at  point  of  stagnation 

underlying  well  , 

(so  -  '. 


V  . 


for  uniform 
sc  reen-ve Iocity 

(1.96) 

1.971 
1..  992 

.  .2.030 
2.087 
2.170 
2.287 
2.459 

« 

r  •  2.733 

;  .  3.264 

!  .  Oo 


for  uniform 
d  raved o vm  at  screen 


2.34 

2.25 


2.15 

(2.25-)a 

2.37 


2.80 

3.70 

oo 


•  ;  ^-n  .present  problem  4m  =  400  feet  1+-  ,,  ’ 

taken  from  line  IV  in  figure  3  is.  0.40  ft  ^  thaJ .  dlstanc®  the  drawdown, 

assumptions  is  chosen  one  gets’either  of  the  f  vT  °n  v/hlch  of  the 
down  at  the.  point  of  stagnation:  ;  ollovang  values  for  the  draw- 

assuming  uniform  screen-velo.city _  .  ...  . 

s0  -  2.09(Q/2rrT)_+  0.40  ft  =(2.09  x'l.94  ft)  *  0.40  ft  =  4.45  ft. 

assuming  uniform  drawdown  at  screen _  V  *  . . 

V  -  2.25(Q/2ffr)  ♦  0.40  ft  -  (2.25  x  1.94  ft)  .  0.40  ft  »  4.76  ft.' 

the  fomBraSdrhar™en1|dopiedafor  y™ 

is  above  the*  4-1,  ^  '-ure  ,-) •  ^OTfcVer,  if  the  pump  intake 

vmte=  insidi  tL  ^  +S°re?n*  °V'ins  to  the  fiction  of  the  upward  moving 

aorJeft-and  W  h"  tL' ^  greatest  at  the  top  of  the 

actlailv  be  nor  e01"’  Ao“rfingly , .-the  screen-velocity  may 

otter  words  aotnaf  ^  than  “  “  Wi'ei1  ^e  drawdown  is  uniform;  In 

set  bv  thesl  *  conditions  probably  lie  somewhere  between  the  two  limits 

sq-c  oy  these  assumptions. 

.  \ 

is  indeDondentnofU+bbedhth6  drav'do'"n  at  the  P°i"t  of  stagnation  theoretically 
that  ratio  d  e  ra  10  of  horizontal  to  vertical  permeabilitv,  provided 

the  sa"a  °  root  a?*iT“.-nfiaite*  &dial  dista“«  are  to  be  divided  by 

.th.  r  ^  ratl°  bef0re  makinS  in  the  drawdown  at 

Other  points.  (In  the  present  problem  that  ratio  is  taken  as  unity.) 
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Figure  4  gives,  the  approximate  drawdown  distribution  along  the  top  and 
bottom  of  the  bed  and  throughout  the  bed  for*40  percent  penetration.  The 
equipotentials  (lines  of  equal  drawdown)  in  that  figure  are  to  be  regarded  as 
diagrammatic,  having  been  sketched  in  by  hand  from  known  values  of  drawdown 
along  the-  top-  and.  bottom  of  the  bed.  ^ 


It  may  bo  seen  that  by  averaging  drawdowns  observed  at  the  top  and  bottom 
of  the  bed  the  approximate  drawdown  for  complete  penetration  is  obtained.  The 
data  from  the  foregoing  example  bear  this  out:  .  .  ...... 


Well 

Observed  s 


top; 


ft 


Computed  sbob,  ft 
Final  adjusted  s,  ft 
Avg.  stop  and  sbot,  ft 


50  percent  penetration  the  drawdowns 


1  v 

5-.02 

5.88 

4.42 

4.45 


2.  .  . 
3.18  . 
3.00 
3  ,09.  . 
■3.09 


3 

2.30 

2.27 

2.28 
2.28+ 


The  difference  is  well  within  the  error  of  plotting  in  this  case.  For 


average  out .exactly* 


It  is  recommended  that  in  running ■  pumping  tests.tq  determine  T  with 
partially  penetrating  wells  the  observation  wells  be. placed  in  pairs,  one  of 
each  pair  tapping  the  top  of  the  bed^  the  other  the  bottom. 


A  remark  about  Slichter’ s  formula  for  a  partially  penetrating  well  may 
not  be  out  of  place  here.  Reference  is  made-  to  Wenzel’ s  equation  (123) 
Water-Supply  Paper  887,  p.  109).  That  formula  is.  unsound  because  it  assumes 
radial  flow  superimposed  upon  hemisph'e'ri-cal-  -flow,,  .it  being  impossible  to 
adjust  the  potential  distribution  satisfa-ot-ori-ly  along  the.ir  common  boundary. 
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